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1.  Introduction 


1.        Introduction 

An  atlas  of  the  Geophysical  Fluid  Dynamics  Laboratory  (GFDL)  Rhomboi- 
dal  30  with  14  Sigma  Levels  (R30S14)  general  circulation  model  (GCM)  is  made 
available  in  order  to  facilitate  the  intercomparison  of  the  climate  and  variability  of 
the  model  with  observations  and  other  modeling  studies.  Our  goal  is  to  reveal  some 
of  the  model's  strengths  and  weaknesses,  providing  a  benchmark  for  future  experi- 
ments with  the  GFDL  model. 

The  R30S14  version  of  the  GFDL  GCM  is  a  global,  spectral,  primitive 
equation  model.  There  are  14  unequally  spaced  sigma  levels  in  the  vertical  (from 
0.9967  to  0.015)  and  a  rhomboidal  truncation  at  wave  number  30,  yielding  a  hori- 
zontal grid  spacing  of  approximately  3.75°  of  longitude  and  2.225°  of  latitude.  The 
model  has  seasonally  varying  insolation,  sea  surface  temperature  (SST)  and  sea  ice, 
but  the  values  are  fixed  for  each  day.  The  SSTs  and  sea  ice  vary  according  to  long- 
term  observed  climatologies  and  the  same  cycle  is  repeated  for  each  year  in  this  17 
year  integration.  The  spectral  representation  of  the  orography  has  been  improved 
by  smoothing  out  some  of  the  artificial  ripples  in  the  field  that  are  generated  by  the 
spectral  transformation.  This  model  also  features  gravity  wave  drag  and  predicted 
clouds.  Stratiform  clouds  form  and  large  scale  precipitation  begins  when  the  rela- 
tive humidity  exceeds  100%.  Subgrid  scale  precipitation  is  parameterized  by  moist 
convective  adjustment.  The  surface  temperature  over  land  is  calculated  assuming 
there  is  no  heat  storage  in  the  ground.  Soil  moisture  is  predicted  using  the  bucket 
method,  in  which  the  ground  can  absorb  up  to  15  cm  of  rainfall  before  runoff  begins 
(Manabe,  1969).  Standard  bulk  aerodynamic  formula  are  used  to  calculate  the  sur- 
face wind  stress  and  sensible  and  latent  heat  fluxes  using  a  constant  transfer  coeffi- 

cient  of  1x10  over  the  oceans  and  3x10  over  land.  More  details  on  the  GFDL 
GCM  can  be  found  in  Gordon  and  Stern  (1982),  Manabe  and  Hahn  (1981)  and  Lau 
(1981). 


2.    Data  Analysis 

The  seasonal  mean  statisties  for  Deeember-January-February  (DJF), 
March- April-May  iMAM).  June-July-August  (JJA),  and  September-October-No- 
\ ember  (SON),  are  calculated  from  16  years  of  daily  output  from  the  GFDL 
R30S  14  model  run.  Data  from  the  model's  sigma  surfaces  are  interpolated  to  pres- 
sure levels:  the  temperature  and  wind  fields  are  interpolated  linearly,  while  specific 
humidity  is  interpolated  logarithmically  to  pressure  surfaces.  In  the  model's  cloud 
cover  scheme,  there  is  no  fractional  cloudiness;  a  grid  point  at  a  given  sigma  level 
is  either  clear  or  completely  overcast.  The  fractional  cloudiness  shown  in  the  atlas 
is  a  measure  of  the  number  of  days  in  the  season  where  cloud  cover  is  present.  For 
the  layer  average  fractional  cloudiness,  the  entire  layer,  made  up  of  2-3  levels,  is 
considered  overcast  if  one  of  the  levels  is  overcast.  Diabatic  heating  is  calculated 
from  the  latent  heat  release  associated  with  condensation,  absorbed  radiation,  ver- 
tical diffusion  and  the  convective  adjustment  process;  both  the  clouds  and  diabatic 
heating  are  shown  on  sigma  surfaces.  The  mass  weighted  diabatic  heating  is  ob- 
tained by  summing  the  diabatic  heating  rate  over  the  14  model  layers. 

Variability  in  the  model  is  diagnosed  using  interannual  standard  deviations 


-      2 
li  =  1 


calculated  from  the  equation:  a(x)  =    |£  (x(-X)   /(n-  1)  ,  where  x,  is  the  average 


for  one  season  of  a  particular  year  (i),  x  is  the  long-term  ( 1 6  year)  seasonal  average, 
and  n  is  the  number  of  years.  Empiracal  Orthogonal  Function  (EOF)  analysis 
(Kutzbach,  1979)  is  used  to  objectively  identify  the  dominant  modes  of  500mb 
height  variability  over  the  Northern  Hemisphere.  The  fist  two  EOFs,  calculated  us- 
ing covariances  rather  than  correlations,  arc  presented. 

The  zonal  average  vertical  cross-sections  are  constructed  from  data  interpo- 
lated to  8  pressure  levels  (1000,  850,  700,  500,  300,  250,  200,  and  100  mb),  except 
for  the  vertical  cross-sections  of  cloud  cover,  mass-stream  function    and  diabatic 


heating,  which  are  constructed  from  the  original  14  sigma  levels  in  the  model.  The 
merdional  mass-stream  function,  \\f ,  is  caculated  following  Peixoto  and  Oort 

l 

(1992),  \i/  =   1*271/? cos  0  [v  •  P]  —  ,  where  R  is  the  radius  of  the  earth,  0  the  lati- 

J  8 

o 

tude,  g  is  gravity,  P  is  the  surface  pressure,  and  the  brackets  denote  a  zonal  average. 
A  solution  for  \j/  is  then  found  by  vertically  integrating  this  equation  starting  at  the 
top  of  the  atmosphere,  where  \j/  =  0 .  The  transient  meridional  momentum,  heat 
and  moisture  fluxes  are  derived  from  the  model  output  using  the  following  equa- 
tions: «V=  («-«)•  (v-v)  ,  v'T  =  (v-v)  •  (T-f)  ,  and  v'q'  =  (v-v)  ■  (q-q)  , 
where  x  is  the  daily  value  and  x  is  the  long-term  average  for  the  appropriate  season. 
These  quantities  are  calculated  daily  and  then  the  seasonal  and  zonal  means  are  con- 
structed from  the  daily  values. 

Time  filtering  of  certain  variables  is  performed  using  a  Lanczos  filter  (Jus- 
tice, 1976;  Duchon,  1979),  with  half  power  weights  at  3-10  days  to  isolate  synoptic 
times,  and  20-100  days  to  examine  the  Madden  and  Julian  Oscillation  (Madden  and 
Julian,  1971).  The  3-10  day  filter  is  also  used  to  examine  the  seasonal  cycle  of  the 
North  Pacific  and  North  Atlantic  storm  tracks  following  Nakamura  (1992).  We 
used  121  weights,  which  translated  to  a  data  loss  of  60  days  at  the  beginning  and 
end  of  the  17  year  time  series.  The  filtering  is  performed  on  daily  anomalies  (rela- 
tive to  the  long-term  monthly  mean)  and  the  total  standard  deviation  or  mean  of 
these  filtered  anomalies  are  then  calculated  (except  in  the  case  of  Figs  113  and  1 14, 
where  the  time  series  of  the  time  filtered  anomalies  are  shown) .  Only  model  years 
2-16  are  used  in  the  standard  deviation  calculations,  since  model  years  1  and  17  are 
truncated  by  the  filtering  process.  For  filtered  covariance  statistics,  each  individual 
variable  (i.e.  v')  is  filtered  before  computing  the  covariances.  A  few  variables, 
such  as  the  precipitation  and  vertical  velocity,  are  spatially  smoothed  using  a  nine 


point  filter  in  order  to  emphasize  large  scale  features.  The  filter  weights  the  central 
point,  the  four  adjacent  points,  and  the  four  corner  points  by  a  ratio  of  1:0.5:0.3. 
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4.  Zonal  Mean  Vertical  Cross-Sections 
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5.  Seasonal  Means  on  Pressure  and  Sigma  Surfaces 
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I  i  jure  23:      700mb  Departure  of  Temperature  from  the  Zonal  Mean 
( '•  C),  Northern  Hemisphere  for  DJF  (top)  and  MAM  (bottom) 
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Figure  24:      700mb  Departure  of  Temperature  from  the  Zonal  Mean 
(°  C),  Northern  Hemisphere  for  JJA  (top)  and  SON  (bottom) 
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ire  25:       700mb  Departure  of  Temperature  from  the  Zonal  Mean 
'       Southern  Hemisphere  for  D.JI;  (top)  and  MAM  (bottom) 
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Figure  26:      700mb  Departure  of  Temperature  from  the  Zonal  Mean 
( °  C),  Southern  Hemisphere  for  JJA  (top)  and  SON  (bottom) 


31 


z*500(m) 


P.ll- 


- 


- 


180  E 


z*50o(m) 


MAM 


O^C  -■ '  ^  -  ^r — s     „  ^~^r 

,-'-- 

'  /     ^!rJ^^>^c^~T^^^^~^-^l 

\\    - 

Y? 

^'—'''k%0^Jvir^     **    ' 

i '-''' 

'///f      J\M   '  zXy{^ 

w 

/ 

HjJjJpffi- 1  oovvyy 

^a&7 

fir 

7     f 

80  H 


175 

150 

125 

100 

75 

50 

25 

0 


"OOmb  Departure  ol'Cieopotential  Height  Iron)  the  Zonal  Mean  (m), 
Northern  Hemisphere  lor  I J>.H   flop)  and  MAM  (bottom) 


32 


z*5oo(m) 


JJA 


3* 


180  E 


500 


(m) 


SON 


80  E 


I 


175 
50 
25 

0 


Figure  28:      500mt>  Departure  of  Geopotential  Height  from  the  Zonal  Mean  irrij. 
Northern  Hemisphere  for  JJA  (top)  and  SON  (bottom) 
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500mb  Departure  of  Geopotential  Height  from  the  Zonal  Mean  (m) 
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Figure  32:      200mb  Departure  of  Geopotential  Height  from  the  Zonal  Mean  (m), 
Northern  Hemisphere  for  JJA  (top)  and  SON  (bottom) 
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Figure  34:      200mb  Departure  of  Geopotential  Height  from  the  Zonal  Mean  (mj 
Southern  Hemisphere  for  JJA  (top)  and  SON  (bottom) 
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Figure  35:      200  mb  Zonal  Wind  (m/s)  for  DJF  (top)  and  MAM  (bottom) 
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Figure  36:      200  mb  Zonal  Wind  (m/s)  for  JJA  (top)  and  SON  (bottom) 
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Figure  38:      200  mb  Meridional  Wind  (m/s)  for  JJA  (top)  and  SON  (bottom) 
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Figure  48:      850  mb  Velocity  Potential  (1x10    rrr/s)  for  JJA  (top)  and  SON  (bottom) 
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Figure  49: 


200  mb  Velocitv  Potent  i;ii  .  i  x  !(/"  nr/si  for  I  J.I  I-  (top)  and  MAM  (bottom) 
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Figure  50:      200  mb  Velocity  Potential  ( lx  10  nr/s)  for  JJA  I  topi  and  SON  (bottom) 
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850  mb  Stream  Function  f  i  x  I  0  '  nr/s)  for  DJF  (top)  and  MAM  (bottom) 
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Figure  52:      850  mb  Stream  Function(  lxlO6  m2/s)  for  JJA  (top)  and  SON  (bottom) 
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F-igure  53: 


200  nib  Stream  Function  ( 1x10    m2/sj  lor  DJF  (top)  and  MAM  (bottom) 
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Figure  54:      200  mb  Stream  Function  ( IxlO6  m2/s)  for  JJA  (top)  and  SON  (bottom) 
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850  mb  Specific  Humidity  (g/kg;  for  DJF  (top)  and  MAM  (bottom) 
icld  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  56:      850  mb  Specific  Humidity  (g/kg)  ior  JJA  (top)  and  SON  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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500  nib  Vertical  Velocity  (mb/day)  for  DJF  (top)  and  MAM  (bottom) 
r!n-  fieid  was  smoothed  using  a9poini  iiltcr 'see  data  analysis  section). 
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Figure  58'      500  mb  Vertical  Velocity  (mb/day)  for  JJA  (top)  and  SON  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 


63 


: 


CLD,      (%) 


[>ll; 


60E  120E  180  120W 

Lon (deg) 


60W 


190 


■  75 


60 
45 
30 


CLDlow  (%) 


MAM 


"... 


■ 


180  120W 

Lon  (deg) 


90 
75 

1  ^^ 

45 
30 


.    :- 

■. 
: 


Low  level  Composite  (sigma  levels  0.997,  0.979,  0.935) 
Clouds  [%)  for  DJF  (top)  and  MAM  (bottom) 
:ld  was  smoothed  using  a  9  poinl  filter  (see  data  analysis  section). 
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Figure  60:      Low  level  Composite  (sigma  levels  0.997,  0.979,  0.935 ) 
Clouds  (%)  for  JJA  (top)  and  SON  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section; 


90 
75 
60 
45 
30 


65 


CLDmid{%) 


n.n 


. 


-- . 


Lon  (deg) 


CLDmid  (%) 


MAM 


Lon  (deg) 


Figure  61 :      Mid  level  Composite  (sigma  levels  0.568,  0.46) 

Clouds  ('/<  )  lor  [III   (top)  and  MAM  (bottom) 

This  field  was  smoothed  using  a  9  point  Idler  (see  data  analysis  section). 
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Figure  62:      Mid  level  Composite  (sigma  levels  0.568.  0.46) 

Clouds  (%)  for  JJA  (top)  and  SON  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  63:       High  level  Composite  (sigma  levels  0.353,  0.257) 
Clouds  (%)  for  DJF  Mopi  and  MAM  (bottom) 
:ld  was  smoothed  using  a  9  point  filter  (see  data  analysis  section 
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Figure  64:      High  level  Composite  (sigma  levels  0.355,  0.257) 

Clouds  (%)  for  JJA  (top)  and  SON  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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figure  65:       Mass  Weighted  Vertical  Integral  of 
)iabatic  Heating  (°C  / 'day)   for  DJF  'top;  and  MAM  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  66:      Mass  Weighted  Vertical  Integral  of 

Net  Diabatic  Heating  (°C/day)   for  JJA  (top)  and  SON  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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6.  Seasonal  Mean  Surface  Plots 


Free  (mm/day)  DJF 


'. 


30N 


_ 


20E  180  120W 

Lon  (deg) 


16 

14 
12 
10 

8 
6 

4 


Prec  (mm/day) 


MAM 


60N 


30N 


W) 


30S- 


20E  180  120W 

Lon  (deg) 


16 
14 
1  12 
I  1  0 

8 


6 


4 


Figure  67-      Precipitation  (mm/day)  lor  DJF  (top)  and  MAM  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  68:      Precipitation  ( mm/day j  for  JJA  (top)  and  SON  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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69        Mean  Sea  Level  Pressure  (mb.i  for  DJF  (top)  and  MAM  (bottom ) 
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Figure  70:      Mean  Sea  Level  Pressure  (mb)  for  JJA  (topj  and  SON  (bottom) 
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f'igure  71        Surface  Wind  Speed  ;m/s;  and  Direction  (equal  length  vectors) 
for  DJF  (top;  and  MAM  (bottom) 
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Figure  72:      Surface  Wind  Speed  im/s>  ana  Direction  (equal  length  vectors; 

for  JJA  (top)  and  SON  (bottom) 
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Figure  73: 


Surface  Wind  Stress  (N/m2)  and  Direction  (equal  length  vectors) 
for  DJF  (top)  and  MAM  (bottom) 
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Figure  74:      Surface  Wind  Stress  (N/m2)  and  Direction  (equal  length  vectors) 

for  JJA  (top)  and  SON  (bottom) 
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Figure  75:      Sensible  Heat  Flux  (W/m2)  for  DJF  (top)  and  MAM  (bottom) 

A  positive  flux  indicates  an  energy  loss  from  the  ocean. 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  76:       Sensible  Heat  Flux  (W/m2)  for  JJA  (top)  and  SON  (bottom) 

A  positive  flux  indicates  an  energy  loss  from  the  ocean. 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  77:      Latent  Heat  Flux  (W/m2)  for  DJF  (top)  and  MAM  (bottom) 

A  positive  flux  indicates  an  energy  loss  from  the  ocean 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  78:      Latent  Heat  Flux  (W/m2)  for  JJA  (top)  and  SON  (bottom) 

A  positive  flux  indicates  an  energy  loss  from  the  ocean. 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  79:      Short  Wave  Radiative  Flux  (W/m2)  for  DJF  (top)  and  MAM  (bottom) 

A  positive  flux  indicates  an  energy  gain  to  the  ocean. 

I  his  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  80:      Short  Wave  Radiative  Flux  (W/m2)  for  JJA  (top)  and  SON  (bottom) 

A  positive  flux  indicates  an  energy  gain  to  the  ocean. 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Long  Wave  Radiative  Flux  (W/m2)  for  DJF  Hop)  and  MAM  (bottom) 
A  positive  flux  indicates  an  energy  loss  from  the  ocean. 
field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  82:      Long  Wave  Radiative  Flux  (W/m2)  for  JJA  (top)  and  SON  (bottom) 

A  positive  flux  indicates  an  energy  loss  from  the  ocean. 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Net  Surface  Heal    Mux  (W/m2)  lor  DJF  (lop)  and  MAM  (bottom) 
A  positive  flux  indicates  an  energy  gain  to  the  ocean. 
This  field  was  smoothed  using  a  ()  point  filter  (see  data  analysis  section). 
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Figure  84:      Net  Surface  Heat  Flux  (W/m2)  for  JJA  (top)  and  SON  (bottom) 

A  positive  flux  indicates  an  energy  gain  to  the  ocean. 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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7.  Interannual  Variability: 
Interannual  Standard  Deviations  and  EOF's 
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Figure  86:      Interannual  Standard  Deviation  of  500  mb  Geopotential  Height  (m) 

for  DJF  (top)  and  MAM  (bottom) 
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Figure  86: 


Interannual  Standard  Deviation  of  500  mb  Geopotential  Height  (m) 
for  DJF  (top)  and  MAM  (bottom) 
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Figure  87:      Interannuai  Standard  Deviation  of  500  mb  Geopotential  Height  (m) 

for  JJA  (top)  and  SON  (bottom) 
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Figure  88:      Interannual  Standard  Deviation  of  200  mb  Zonal  Wind  (m/s) 
for  DJF  (top)  and  MAM  (bottom) 
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Mgure  89:      Interannual  Standard  Deviation  of  200  mb  Zonal  Wind  (m/s) 
for  JJA  (top)  and  SON  (bottom) 
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Figure  90:      Interannual  Standard  Deviation  of  Sea  Level  Pressure  (mb) 
for  DJF  (top)  and  MAM  (bottom) 
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Figure  9 1 :      Interannual  Standard  Deviation  of  Sea  Level  Pressure  ( mb  I 
for  JJA  (top;  and  SON  (bottom) 
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Figure  92:      500mb  Height  Empirical  Orthogonal  Functions 
(using  covariances)  for  EOF  1  (top)  and  EOF  2  (bottom) 
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Figure  93:      500mb  Height  Empirical  Orthogonal  Functions 
(using  covariances)  for  EOF  1  (top)  and  EOF  2  (bottom) 
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Figure  94:      500mb  Height  Empirical  Orthogonal  Functions 
(using  covariances)  for  EOF  1  (top)  and  EOF  2  (bottom) 
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Figure  95:      500mb  Height  Empirical  Orthogonal  Functions 
for  EOF  1  (top j  and  EOF  2  (bottom) 
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8.  Middle  Latitude  Processes 
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Figure  96:      Mean  3-10  Day  Filtered  850mb  Transient  Meridional  Heat 

Flux  [°C  ■  m/s)  for  DJF  (top)  and  MAM  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  97:      Mean  3-10  Day  Filtered  850mb  Transient  Meridional  Heat 

Flux  (°C  •  m/s)  for  JJA  (top)  and  SON  (bottom) 
This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  98:       Mean  3-10  Day  Filtered  850mb  Transient  Meridional  Moisture 

Flux  (g  ■  m/kg  ■  s)  for  DJF  (top)  and  MAM  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  99:      Mean  3-10  Day  Filtered  850mb  Transient  Meridional  Moisture 

Flux  (#  •  m/kg  ■  s)  for  JJA  (top)  and  SON  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  100:      Mean  3-10  Day  Filtered  200mb  Transient  Meridional  Flux  of  Zonal 

Momentum  (m2/s2)  for  DJF  (top)  and  MAM  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  101:      Mean  3-10  Day  Filtered  200mb  Transient  Meridional  Flux  of  Zonal 

Momentum  (m2/s2)  for  JJA  (top)  and  SON  (bottom) 

This  field  was  smoothed  using  a  9  point  filter  (see  data  analysis  section). 
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Figure  102:      Mean  3-10  Day  Filtered  200mb  Kinetic  Energy  (m2/s2) 
for  DJF  (top)  and  MAM  (bottom) 
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Figure  103:      Mean  3-10  Day  Filtered  200mb  Kinetic  Energy  (m2/s2) 
for  JJA  (top)  and  SON  (bottom) 
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Figure  1 10:      3-10  Day  Filtered  Standard  Deviation  of  250mb  Geopotential 

Heignt  (m)  (top)  and  Monthly  Mean  250mb  Zonal 

Wind  (m/s)(bottom),  averaged  30N-50N. 
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9.  Tropical  Processes 
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Figure  115:     250mb  Velocity  Potential  (  1x10   m2/s)  filtered  to  retain  periodicities  of 
20-100  days,  averaged  ION- 1  OS,  for  model  year  4 
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Figure  1 16:     250mb  Velocity  Potential  (1x10    m2/s)  filtered  to  retain  periodicities  of 
20-100  days,  averaged  10N-10S.  for  model  year  7 
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